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Abstract 

Consider the Klein-Gordon equation (KGE) in IR", n > 2, with constant or variable 
coefficients. We study the distribution of the random solution at time t € IR. We 
assume that the initial probability measure fiQ has zero mean, a translation-invariant 
covariance, and a finite mean energy density. We also asume that satisfies a Rosenblatt- 
or Ibragimov-Linnik-type mixing condition. The main result is the convergence of fit to 
a Gaussian probability measure as t ^ oo which gives a Central Limit Theorem for the 
KGE. The proof for the case of constant coefficients is based on an analysis of long time 
asymptotics of the solution in the Fourier representation and Bernstein's 'room-corridor' 
argument. The case of variable coefficients is treated by using an 'averaged' version of 
the scattering theory for infinite energy solutions, based on Vainberg's results on local 
energy decay. 
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1 Introduction 



The aim of this paper is to underhne a special role of equilibrium distributions in statistical 
mechanics of systems governed by hyperbolic partial differential equations (for parabolic equa- 
tions see [6, 27]). Important examples arise when one discusses the role a canonical Gibbs 
distribution (CGD) in the Planck theory of spectral density of the black-body emission and 
in the Einstein-Debye quantum theory of sohd state (see, e.g. [31]). [The word 'canonical' is 
used in this paper to emphasize the fact that the probability distribution under consideration is 
formally related to the 'Hamiltonian', or the energy functional, of the corresponding equation 
by the Gibbs exponential formula. Owing to the linearity of our equations, there are plenty 
of other first integrals which lead to other stationary measures.] Historically, the emission 
law was estabhshed at a heuristical level by Kirchhoff in 1859 (see [34]) and stated formally 
by Planck in 1900 (see [25]). The law concerns the correspondence between the temperature 
and the colour of an emitting body (e.g., a burning carbon, or an incandescent wire in an 
electric bulb). Furthermore, it provides a fundamental information on an interaction between 
the Maxwell field and a 'matter'. Planck's formula specifies a 'radiation intensity' It{^) of 
the electromagnetic field at a fixed temperature T > 0, as a function of the frequency a; > 0. 
It is convenient to treat /t(") as the spectral correlation function of a stationary random pro- 
cess. Then if denotes an equilibrium distribution of this process, the Kirchhoff- Planck law 
suggests the long-time convergence 

/^t ^ QT: oo. (1.1) 

Here fit is the distribution at time t of a nonstationary random solution. The resulting equi- 
librium temperature T is determined by an initial distribution /iq. Convergence to equilibrium 
(1.1) is also expected in a system of Maxwell's equations coupled to an equation of evolution 
of a 'matter'. For example, both (1.1) and the Kirchhoff-Planck law should hold for the cou- 
pled Maxwell-Dirac equations [5], or for their second-quantised modifications. However, the 
rigorous proof here is still an open problem. 

Previously, the convergence of type (1.1) to a GGD has been established for an ideal gas 
with infinitely many particles by Sinai (see, e.g., [7]). Similar results were later obtained for 
other infinite-dimensional systems (see [2, 10] and a survey [9]). For nonlinear wave problems, 
the first result of such kind has been estabhshed by Jaksic and Pillet in [18]: they consider 
a system of a classical particle coupled to a wave field in a smooth nonlocal fashion. For all 
these models, the CGD gx is well-defined, although the convergence is highly non-trivial. On 
the other hand, for the local coupling such as in the Maxwell-Dirac equations, the problem 
of 'ultraviolet divergence' arises: the CGDs cannot be defined directly as the local energy 
is formally infinite almost surely. This is a serious technical difficulty that suggests that, to 
begin with, one should analyse convergence to non-canonical stationary measures /loo, with 
finite mean local energy: 

fit fiocj t ^ oo. (1.2) 

In fact, most of the above-mentioned papers establish the convergence to both CGDs and 
non-canonical stationary measures, by using the same methods. In our situation, the afore- 
mentioned ultraviolet divergence makes the difference between (1.1) and (1.2). 

In this paper we prove convergence (1.2) for the Klein-Gordon equation (KGE) in IR", 
n>2: 

= E]=i{dj-iAj{x)fu{x,t)-m^u{x,t), x e IR", 
Uo{x), u\t=o = Vo{x). 
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Here dj = — — , x e IR", i e IR, m > is a fixed constant and {Ai(x), . . . ,An(x)) a vector po- 

OXj 

tential of an external magnetic field; we assume that functions Aj (x) vanish outside a bounded 
domain. The solution u{x,t) is considered as a complex-valued c/asszca/ function. 

It is important to identify a natural property of the initial measure /xq guaranteeing conver- 
gence (1.2). We follow an idea of Dobrushin and Suhov [10] and use a 'space'-mixing condition 
of Rosenblatt- or Ibragimov-Linnik-type. Such a condition is natural from physical point of 
view. It replaces a 'quasiergodic hypothesis' and allows us to avoid introducing a 'thermostat' 
with a prescribed time-behaviour. Similar conditions have been used in [2, 3, 33, 32]. In this 
paper, mixing is defined and applied in the context of the KGE. 

Thus we prove convergence (1.2) for a class of initial measures /xq on a classical function 
space, with a finite mean local energy and satisfying a mixing condition. The limiting measure 
/ioo is stationary and turns out to be a Gaussian probability measure (GPM). Hence, this result 
is a form of the Central Limit Theorem for the KGE. 

Another important question we discuss below is the relation of the limiting measure fi^o to 
the CGD qt- The (formal) Klein-Gordon Hamiltonian is given by a quadratic form and so the 
CGDs qt are also GPMs, albeit generalised (i.e. hving in generahsed function spaces). As our 
limiting measures fi^o are 'classical' GPMs. they do not include CGDs. However, in the case 
of constant coefficients, a CGD can be obtained as a limit of measures as the 'correlation 
radius' figuring in the mixing conditions imposed on /^q tends to zero. More precisely, we 
assume that for a fixed T > 

^E(vo{x)vo{y) + Vuoix)-Vuo{y)+m\oix)uoiy))^T5{x-y), r ^ 0, (1.4) 

where E denotes the expectation. Then the covariance functions (CFs) of corresponding limit 
GPM converge to the covariance functions of the CGD qt- In turn, this implies the 
convergence 

A*t ^ A«oo ~ 9t, r^l. (1.5) 

See Section 4. 

It should be noted that the existence of a 'massive' (in a sense, infinite-dimensional) set of 
the limiting measures fi^o that are different from CCD's is related to the fact that KGE (1.3) is 
degenerate and admits infinitely many 'additive' first integrals. Like the Klein-Gordon Hamil- 
tonian, these integrals are quadratic forms; hence they generate GPMs via Gibbs exponential 
formulas. 

Convergence (1.2) has been obtained in [19, 20, 21] for translation-invariant initial mea- 
sures ^Q. However, the original proofs were too long and used a specific apparatus of Bessel's 
functions applicable exclusively in the case of the KGE. They have not been published in detail 
because of lack of a unifying argument that could show the limits of the method and its forth- 
coming developments. To clarify the mechanism behind the results, one needed some new and 
robust ideas. The current work provides a modern approach applicable to a wide class of linear 
hyperbolic equations with a nondegenerate 'dispersion relation', sec Eqn (7.20) below. We also 
weaken considerably the mixing condition on measure /iq. Moreover, our approach yields much 
shorter proofs and is applicable to non-translation invariant initial measures. The last fact is 
important in relation to two-temperature problem [3, 12, 33] and hydrodynamic limit [8]. Such 
progress became possible in a large part owing to a systematic use of a Fourier transform 
(FT) and a duality argument of Lemma 7.1. [The importance of the Fourier transform were 
demonstrated in earlier works [3, 32, 33].] 
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Similar results, for the wave equation (WE) in IR" with odd n > 3, are established in [11] 
which develops the results [26]. The KGE shares some common features with the WE (which 
is formally obtained by setting m = in (1.3)), and the exposition in [20, 21] followed the 
structure of the earher work [26]. On the other hand, the KGE and WE also have serious 
differences, see below. 

It is worth mentioning that possible extensions of our methods include, on the one hand, 
Dirac's and other relativistic-invariant linear hyperbolic equations and on the other hand har- 
monic lattices, as well as 'coupled' systems of both types. We intend to return to this problems 
elsewhere. 

We now pass to a detailed description of results. Formal definitions and statements are 
given in Section 2. Set: Y{t) = {Y%t),Y^{t)) = {u{- , t) , u{- , t)) , Fq = {Yo^,Y^') = {uo,Vo). 
Then (1.3) takes the form of an evolution equation 

Y(t) = AY(t), t e IR; y(0) = Yo. (1.6) 

Here, 

A- (I I), 



where A — jyj=i{dj — iAj{x)y — m^. We assume that the initial date Yq is a random element of 
a complex functional space Ti corresponding to states with a finite local energy, see Definition 
2.1 below. The distribution of Yq is a probability measure of mean zero satisfying some 
additional assumptions, see Conditions S1-S3 below. Given t G IR, denote by iit the measure 
that gives the distribution of Y{t), the random solution to (1.6). We study the asymptotics of 
/It SiS t ^ ±oo. 

We identify (D = IR^ and denote by tensor product of real vectors. The CPs of the initial 
measure are supposed to be translation- invariant: 

Qi^{x,y):^E{Y^{x)^Y^{y))^qi^{x-y), x,yeJR\ i,j = 0,l, (1.8) 

(in fact our methods require a weaker assumption, but to simplify the exposition, we will not 
discuss it here). We also assume that the initial mean energy density is finite: 

eo := E{\vo{x)\'' + \Vuo{x)\^+m^\uo{x)\^) = g^'(O)- Ag°°(0) ^^^(O) < oo, x e IR". (1.9) 

Finally, we assume that measure /lo satisfies a mixing condition of a Rosenblatt- or Ibragimov- 
Linnik type, which means that 

Yq{x) and lo(y) are asymptotically independent as \x — y\ oo. (1-10) 

As was said before, our main result gives the (weak) convergence (1.2) of fit to a limiting 
measure jj,^ which is a stationary GPM on 7i. A similar convergence holds for t — > — oo. An 
explicit formulas are then given for the CFs of /Xqo- 

The strategy of the proof is as follows. First, we prove (1.2) for the equation with constant 
coefficients {Ak{x) = 0), in three steps. 

I. We check that the family of measures /z^, i > 0, is weakly compact. 

II. We check that the CFs converge to a hmit: for i,j — 0, 1, 

Qnx,y)= jY\x)®Y^{y)fit{dY)^Q%{x,y), t ^ oo. (1.11) 
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III. Finally, we check that the characteristic functionals converge to a Gaussian one: 



M^) := J exp{t{Y, ^)}MdY) ^ exp{-^Qoo(^, ^)}, t ^ oo. (1.12) 

Here \& is an arbitrary element of the dual space and Qoo the quadratic form with the integral 
kernel {Q^Joi^^ y))i,j=o,i'-, 0^^ ^) denotes scalar product in a real Hilbert space L^(1R"') (8) IR^. 

Property I follows from the Prokhorov Theorem by a method used in [37]. First, we prove a 
uniform bound for the mean local energy in fit, using the conservation of mean energy density. 
The conditions of the Prokhorov Theorem are then checked by using Sobolev's embedding 
Theorem in conjunction with Chebyshev's inequality. Next, we deduce property II from an 
analysis of oscillatory integrals arising in the FT. An important role is attributed to Proposition 
6.1 reflecting the properties of the CFs in the FT deduced from the mixing condition. 

On the other hand, the FT approach alone is not sufficient for proving property III even in 
the case of constant coefficients. The reason is that a function of infinite energy corresponds 
to a singular generalised function in the FT, and the exact interpretation of mixing condition 
(1.10) for such generalised functions is unclear. We deduce property III from a representation 
of the solution in terms of the initial date in the coordinate space. This is a modification of 
the approach adopted in [19, 20, 21]. It allows us to combine the mixing condition with the 
fact that waves in the coordinate space disperse to infinity. This leads to a representation 
of the solution as a sum of weakly dependent random variables. Then (1.12) follows from a 
Central Limit Theorem (CLT) under a Lindeberg-type condition. Checking such a condition 
is an important part of the proof. 

It is useful to discuss the dispersive mechanism that is behind (1.12) and compare the KGE 
(m > 0) and WE (m = 0). Take, for simplicity, n = 3 and Mq = 0. The solution to (1.3) (with 
Ak{x) = 0) is given by 

u{x,t)= £{x-y,t)vo{y)dy, t > 0, (1-13) 



where S is the 'retarded' fundamental solution 

Ji is the Bessel function of the first order. For m — function E{-.,t) is supported by the 
sphere \x\ — t oi area ~ t"^, and (1.13) becomes the Kirchhoff formula 

u{x,t)^^^ J vo{y)dS{y), (1.15) 

\x-y\=t 

which manifests the dispersion of waves in the 3D space. Dividing the sphere {y e IR^ : \x — 
y\ — t} into N t'^ 'rooms' of a fixed width d » 1, we rewrite (1.15) as 

N 

u{x,t)^^, (1.16) 



where are nearly independent owing to the mixing condition. Then (1.2) follows by well- 
known Bernstein's 'room-corridor' arguments. 
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For m > function S{-,t) is supported by the ball < t which means the absence of a 
strong Huyghen's principle for the KGE. The volume of the ball is ~ t^, hence rewriting (1.13) 
in the form (1.16) would need asymptotics of the type 

S{x,t) = 0{t-^/^), \x\<t (1.17) 

as i — > oo. As Ji(r) ~ cos(r — 37r/4)/y^, asymptotics (1-17) only holds in the region \x\ < vt 
with V < 1. For instance, 

I cos(m7t — 37r/4) 

^\\x\=vt ' (^i)3/2 ' 



where 7 = \/l — v'^. However, the degree of the decay is different near the light cone = t 
corresponding to v = 1 and 7 = 0. For example, for a fixed r > 0, 



p\ cos{my/2H - 37r/4) 3/4^ ^o^ 



where r — t — \x\ is the 'distance' from the hght cone. This illustrates that an apphcation of 
Bernstein's method in the case of the KGE requires a new idea. 

The key observation is that the asymptotics (1.18) displays oscillations ~ cosm\/2rt of £ 
near light cone as t — > oo. The solution becomes an oscillatory integral, and one is able to 
compensate the weak decay ~ t~^/'^ by a partial integration with Bessel functions, by method 
following an argument from [23, Appendix B]. Such an approach was used in [21] and was 
accompanied by tedious computations in a combined 'coordinate- momentum' representation. 
The approach adopted in this paper allows us to avoid this part of the argument. An important 
role plays a duality argument of Lemma 7.1 leading to an analysis of an oscillatory integral 
with a phase function (='dispersion relation') with a nondegenerate Hessian, see (7.20). 

Simple examples show that the convergence may fail when the mixing condition does not 
hold. For instance, take Uq{x) = ±1 and Vo{x) = with probability p± = 0.5 then mean value 
is zero and (1.9) holds, but (1.10) does not. The solution u{x, t) = ±cos {mt) a.s., hence /it 
is periodic in time, and (1.2) fails. 

Finally, a comment on the case of variable coefficients Ak{x). In this case explicit formulas 
for the solution are unavailable. Here we construct a scattering theory for solutions of infinite 
global energy. This version of the scattering theory allows us to reduce the proof of (1.2) to 
the case of constant coefficients (this strategy is similar to [4, 11, 12]). In particular, in [11] 
one establishes, in the case of a WE, a long-time asymptotics 

U{t)Yo = eUo{t)Yo + p{t)Yo, t>0. (1.19) 

Here U{t) is the dynamical group of the WE with variable coefficients, Uo{t) corresponds to 
'free' equation, with constant coefficients, and G is a 'scattering operator'. In this paper, 
instead of (1.19), we use a dual representation: 

U'(t)^ ^Uo(t)W^ + r(t)'^, t>0. (1.20) 

Here U'{t) is a 'formal adjoint' to the dynamical group of Eqn (1.3), while UQ{t) corresponds 
to the 'free' equation, with Ak{x) = 0. The remainder r{t) is small in mean: 

£;|(ro,r(t)*)|'^0, t^oo. (1.21) 
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This version of scattering theory is essentially based on Vainberg's bounds for the local energy 
decay (see [35, 36]). 

Remark i) In [11] we deduce asymptotics (1.19) from its primal counterpart (1.20). In this 
paper we do not analyse connections between (1.20) and (1.19). 

ii) It is useful to comment on the difference between two versions of scattering theory produced 
for the WE and KGE. In the first theory, the remainders p{t) and r{t) are small a. s., while in 
the second theory, developed in this paper, r(t) is small in mean (see (1.21)). Such a difference 
is related to a slow (power) decay of solutions to the KGE. 

The main result of the paper is stated in Section 2 (see Theorem A). Sections 3-8 deal with 
the case of constant coefficients: the main statement is given in Section 3 (see Theorem B), 
the relation to CGDs is discussed in Section 4, the compactness (Property I) is established in 
Section 5. convergence (1.11) in Section 6, and convergence (1.12) in Sections 7, 8. In Section 
9 we check the Lindeberg condition needed for convergence to a Gaussian limit. In Section 10 
we discuss the infinite energy version of the scattering theory, and in Section 11 convergence 
(1.2). In Appendix A we collected FT-type calculations. Appendix B concerns with a formula 
on generalised GPMs on Sobolev's spaces. 

Acknowledgements The authors thank V.I. Arnold, A. Bensoussan, LA. Ibragimov, 
H.P. McKean, J. Lebowitz, A.I. Shnirelman, H. Spohn, B.R. Vainberg and M.I. Vishik for 
fruitful discussions and remarks. 



2 Main results 
2.1 Notation 

We assume that functions Ak{x) in (1.3) satisfy the following conditions: 

El. Aj[x) are real C°^-functions. 

E2. Aj{x) = for \x\ > Rq, where Rq < oo. 

0X2 OX I 

Assume that the initial date Iq belongs to the phase space Ti defined below. 

Definition 2.1 H = i^/„^(IR") © i7i°^(lR") is the Frechet space of pairs Y{x) = {u{x),v{x)) 
of complex functions u{x), v{x), endowed with local energy seminorms 

\\Y\\\= J [\v{x)\'^ + \Vu{x)\'^ + m'^\u{x)\'^)dx < 00, Vi? > 0. (2.1) 

|a;|<R 

Proposition 2.2 follows from [22, Thms V.3.1, V.3.2]) as the speed of propagation for Eqn 
(1.3) is finite. 

Proposition 2.2 i) For any Yq & 7i there exists a unique (generalised) solution Y{t) e 
C{M,n) to (1.6). 

ii) For any t e IR the operator U{t) : Yq ^ Y{t) is continuous in Ti. 
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Let us choose a function ({x) G Co~(]R") with C(O) 7^ 0. Denote by i/f„^(]R"), s e M, the 
local Sobolev spaces, i.e. the Prechet spaces of distributions u e D'(1R") with finite seminorms 

\\u\\s,R := ||A^(C(x/i?)xi)|U2(Hn), (2.2) 

where A''^v := Ff7^^{{k)^v{k)), (k) := \J\k\'^ + 1, and v := Fv is the FT of a tempered distri- 
bution V. For -0 e -D define Fi/j{k) = e^^''''il){x)dx. 

Definition 2.3 For s G IR denote = Hl^^^lBJ") © i/foc(lR")- 

Using standard techniques of pseudodifferential operators (see, e.g. [16]) and Sobolev's 
Theorem, it is possible to prove that TiP — H C TiT^ for every £ > 0, and the embedding is 
compact. 

2.2 Random solution. Convergence to equilibrium 

Let (f2, S, P) be a probability space with expectation E and i3(7Y) denote the Borel o"-algebra 
in Ti. We assume that Yq = Yq{uj,-) in (1.6) is a measurable random function with values 
in {7i, B{7i)). In other words, (a;, a;) 1— >• Yo{u,x) is a measurable map Q x IR" — >• (D^ with 
respect to the (completed) cr-algebras E x i3(lR") and i3((D^). Then, owing to Proposition 2.2, 
Y{t) = U{t)Yo is again a measurable random function with values in {H,B{'H)). We denote 
by ^oldYo) a probability measure on 7i giving the distribution of the Yq. Without loss of 
generality, we assume (il, E, P) = (TC, B{T-C), Ho) and Yq{u!,x) = u!{x) for ^Q^duj) x dx-almost 
all (a;, a;) G 7f x IR". 

Definition 2.4 is a probability measure on Ti. which gives the distribution ofY{t): 

lit{B)= iio{U{-t)B), BeB{n), telR. (2.3) 

Our main goal is to derive the weak convergence of the measures /it in the Frechet space 
for each £ > 0, 

l^t — ^ /^oo, 00, (2.4) 

where //oo is a limiting measure on the space Ti.. This means the convergence 

j f{Y)fXt{dY) ^ J f{Y)fx^{dY), t ^ 00, (2.5) 

for any bounded continuous functional / on 7Y^^. Recall that we identify (D = IR^ and ® 
stands for tensor product of real vectors. Denote = IR^ IR^. 

Definition 2.5 The CPs of the measure /it are defined by 

Q\\x, y) = e{y\x, t) ® Y^{y, t)), i,j = 0, 1, for almost all x,yeWCx IR", (2.6) 

assuming that the expectations in the RHS are finite. 
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WesetP = D®D, and (F, ^) = (Y^,^>^) + {Y^,^>^) for F = G 7^ and ^ = e 

X>. For a probability measure ji on denote by fx the characteristic functional (the FT) 



A probability measure ii is called a GPM (of mean zero) if its characteristic functional has the 
form 

/i(^) = exp{-iQ(^,^)}, *eP, 

where Q is a real nonnegative quadratic form in P. A measure // is called translation- invariant 
if 

liiThB) = /i(5), B G B{H), h e IR", 
where ThY{x) = Y{x - h), x e M". 

2.3 Mixing condition 

Let 0{r) denote the set of all pairs of open bounded subsets A,Bg 1R" at distance dist(^, B) > 
r and (t{A) the a- algebra in Ti generated by the linear functionals F i— > {Y,"^) where \I' G 7? 
with supp ^ C .4. Define the Ibragimov-Linnik mixing coefficient of a probability measure /iq 
on n by (cf [17, Dfn 17.2.2]) 

ip{r) = sup sup — . (2.7) 

(AB)eO(r) ^ g 5 ^ ^(g) MB) 

l^o{B) > 

Definition 2.6 T/ie measure /lo satisfies the strong, uniform Ibragimov-Linnik mixing condi- 
tion if 

(p{r) ^0, r ^ oo. (2.8) 
Below, we specify the rate of decay of (p (see Condition S3). 

2.4 Main assumptions and results 

We assume that measure no has the following properties S0-S3: 

50 Ho has zero expectation value, 

EYo{x) = 0, xe IR". (2.9) 

51 iiQ has translation-invariant CFs. i.e. Eqn (1.8) holds for almost all x,y E IR". 
32/10 has a finite mean energy density, i.e. Eqn (1.9) holds. 

S3 //o satisfies the strong uniform Ibragimov-Linnik mixing condition with 

oo 

lp = J r''-^ip^/'^(r)dr < oo. (2.10) 
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Define, for almost all x,y ^ IR"^, the matrix Qoo{x, y) = (Q'^J^{x, yj^ . . ^ ^ by 

O (rv)^U ^< + ^ * - (^o' - - (2 11) 

Qo.{x,y) _ 2 (^10 _ ^oi)(^ _ (^11 _ _ m')qOO^{x - y) ) ' ^''^^^ 

Here V{z) is the fundamental solution for the operator —A + m^, and * stands for the con- 
volution of generahzed functions. We show below that q^^ e L^(IR") (see (6.1)). Then the 
convolution V * q^^ in (2.11) also belongs to L^(IR"'). 

Let H = L^(1R") © i?^(IR") denote the space of complex valued functions * = (*o, *i) 
with a finite norm 

11*11^ = / (l*o(2;)|' + |V*i(x)|' + |*i(x)|') dx < oo. (2.12) 

Denote by Qoo a real quadratic form in H defined by 

Qoo{^,^)= E J {Q'iix,y)^,{x),^,{y))dxdy, (2.13) 

where (^-j stands for real scalar product in (D^ = IR''. The form is continuous in H by 
Corollary 6.3. 

Theorem A Let n > 2, m > 0, and assume that El— E3, SO— S3 hold. Then 

i) The convergence in (2.4) holds for any £ > 0. 

ii) The limiting measure /i^ is a GPM on H. 

Hi) The characteristic functional of /i^ has the form 

where W -.V ^ H is a linear continuous operator. 



2.5 Remarks on conditions on the initial measure 

i) The (rather strong) form of mixing in Definition 2.6 is motivated by two facts: (a) it greatly 
simplifies the forthcoming arguments, (b) it allows us to produce an 'optimal' (most slow) 
decay of 99 indicating natural limits of Bernstein's room-corridor method. Condition (2.7) can 
be easily verified for GPMs with finite-range dependence and their images under 'local' maps 
H ^ H. See the examples in Section 2.6 below. 

ii) The uniform Rosenblatt mixing condition [30] also suffices, together with a higher power 
> 2 in the bound (1.9): there exists 6 > such that 

E^lvoix)]^"-' +\Vuo{x)\'^' + m'\uoix)\'+') < 00. (1.4') 

Then (2.10) requires a modification: 

roo ^ \ 

/ r"'^^a^{r-)dr < 00, where p = min(- -, -), (2.10') 

Jo 2 -|- 2 

where a{r) is the Rosenblatt mixing coefficient defined as in (2.7) but without //(-B) in the 
denominator. The statements of Theorem A and their proofs remain essentially unchanged, 
only Lemma 8.2 requires a suitable modification [17]. 
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2.6 Examples of initial measures with mixing condition 



2.6.1 Gaussian measures 

In this section we construct initial GPMs /iq satisfying SO — S3. Let /iq be a GPM in the 
space H with the characteristic functional 

/ioW=^exp(z(r,^)) = exp{-^Qo(^,^)}, "^eV. (2.14) 

Here Qo is a real nonnegative quadtratic form with an integral kernel {Q o {x,y))i j=q^i. Let 

Q'^{x,y)^qi,\x-y), (2.15) 

for any i, j, where the function Qq e C^(1R") has a compact support. Then SO, SI and 
S2 are satisfied; S3 holds with ip{r) = for r > ro if qo{z) = for j^j > tq. For a given 
matrix function (^q^{zyj such a measure exists in the space Ti iff the corresponding FT is a 

nonnegative matrix-valued measure: (qo{k)^ >0,ke IR", [15, Thm V.5.1]. For example, all 
these conditions hold if qo{k) — Di5'^^f{ki) • . . . • f{kn) with -Dj > and 



2.6.2 Non- Gaussian meaisures 



Now choose a pair of odd functions G C^(IR), with bounded first derivatives. Define //q 

as the distribution of the random function (/°(y°(x)), p{Y^{x))), where (y°, Y^) is a random 
function with a Gaussian distribution from the previous example. Then S0-S3 hold for fi^ 
with a mixing coefficient ^p*{r) = for r > tq. Measure /Xq is not Gaussian if > and the 
functions /* are bounded and nonconstant. 



3 Equations with constant coefficients 

In Sections 3-9 we assume that coefficients Ak{x) = 0. Problem (1.3) then becomes 

{u{x, t) = Au{x, t) — rn^u^x, t), t e IR, 
u\t=o = Uo{x), u\t=o = Vo{x). 



(3.1) 



As in (1.6), we rewrite (3.1) in the form 

Y{t)=AoY{t), telR; Y{0) = Yo. (3.2) 

Here we denote 

Ao-('ll). (3.3) 

where Aq = A — w?. Denote by f/o(t), t G IR, the dynamical group for problem (3.2), then 
Y{t) — UQ{t)YQ. The following proposition is well-known and is proved by a standard partial 
integration. 
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Proposition 3.1 Let Yq = {uo,vo) G H, and Y{-,t) = {u{-,t),u{-,t)) e C{1R,H) is the 
solution to (3.1). Then the following energy hound holds: for R> and t 

J [\u{x,t)\'^ + \Vu(x,t)\^+m'^\u{x,t)\^')dx < J [\vo(x)^\ + \Vuo(x)\'^+m'^\uo(x)\^)dx.(3.A) 

\x\<R \x\<R+\t\ 

Set i^tiB) = no{Uo{-t)B), B e t e IR. Then our main result for problem (3.2) is 

Theorem B Let n>l,m>Q, and Conditions SO— S3 hold. Then the conclusions of Theorem 
A hold with W — I , and limiting measure /loo is translation-invariant. 

Theorem B can be deduced from Propositions 3.2 and 3.3 below, by the same arguments as in 
[37, Thm XII.5.2]. 

Proposition 3.2 The family of measures {iJ,t,t G IR}, is weakly compact in 7i~'^ with any 
£ > 0, and the bounds hold: 

sup E\\Uo{t)Yo\\R < oo, R>0. (3.5) 
Proposition 3.3 For every ^ eV, 

fiti^) = J exp(z(F, ^))i^t{dY) ^ exp{-iQoo(^, ^)}. t ^ oo. (3.6) 

Propositions 3.2 and 3.3 are proved in Sections 5 and 7-9, respectively. We will use repeatedly 
the FT (12.2) and (12.3) from Appendix A. 



4 Relation to CCDs 

In this section we discuss how our results are related to CGDs. We restrict consideration to 
the case of Eqn (1.3) with constant coefficients and to the translation-invariant isotropic case. 
The CGD gx with the absolute temperature T > is defined formally by 

H 

gridu X dv) = — e T Y\_du{x)dv{x), (4.1) 

where H := J (W.)|^ + |V„(.)r + ,n>(.)|>x. and Z is a normalisation eonstant. 

To make the definition rigorous, let us introduce a scale of weighted Sobolev's spaces 
i7*'"(lR'^) with arbitrary s, a G IR. We use notation (2.2). 

Definition 4.1 i) i7*'°(IR"') is the complex Hilbert space of the distributions w e 5"(IR") with 
the finite norm 

\\w\\s,a = ||(x)°A'w||L2(Rn) < oo. (4.2) 
a) H^'°' is the Hilbert space of the pairs Y — {u,v) e i?^"'"*'"(IR") © H^'°'{1R"') with the norm 

\\\Y\\\s,a = \\u\\l+s,a + \\v\\s,a- (4.3) 
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Note that C Ti*'" if s < s and a < a, and this embedding is compact. These facts 
follow by standard methods of pseudodifferential operators and Sobolev's Theorem (see, e.g. 
[16]). 

Now we can define the CGDs rigorously: qt is a GPM on a space W'"', s,a < —n/2, with 
the CFs 

g^%x -y)^ TV{x - y), g].\x - y) ^ T5{x - y), g^\x - y) ^ g}p{x - = 0. (4.4) 

By Minlos Theorem [15, Thm V.5.1], such a measure exists on with s, a < —n/2 as, 
formally (see Appendix B), 

' \Yll^ gridY) < oo. (4.5) 



Measure gx is stationary for the KGE, as its CFs are stationary; the last fact follows from 
formulas (12.6), (12.2). Also, gx is translation invariant, so SI holds. Condition S2 fails since 
the 'mean energy density' fif^ (0) — A5f^°(0) + m'^g^{0) is infinite; this gives an 'ultraviolet 
divergence'. Mixing condition S3 holds due to an exponential decay of the V{z). The conver- 
gence of type (1.1) holds for initial measures /xq that are absolutely continuous with respect 
to the CGD gT, and the limit measure coincides with gx- This mixing property (and even 
if-property) can be proved by using well-known methods developed for Gaussian processes [7], 
and we do not discuss it here. 

Remark Assumption S2 implies that iioiTi) = 1 and hence jiooO^) = 1- This excludes the 
case of a limiting CGD as it is a generalised GPM not supported by H. However, it is possible 
to extend our results to a class of generalised initial measures converging to CGDs. For the 
case of constant coefficients such an extension could be done by smearing the initial generalised 
field as the dynamics commutes with the averaging (cf. [12]). For variable coefficients such 
extension requires a further work. 

To demonstrate the special role of the CGDs we consider a family of initial GPMs ^o,r, 
r e (0, 1], satisfying S0-S3, with the radius of correlation r. More precisely, suppose that the 
corresponding CFs Qq .^ have the following properties G0-G3: 

GO <(^)=<(-^), zeJR\ 

Gl q'^^iz) - Agoo (^) + m'q^^iz) = 0, \z\ > r. 

G2 For some T > 0, ^ J [ ql^,{z) - Ag°° (z) + m\^^,{z) )dz^T, r ^ 0. 

sup / (|<U^)l + |A<°(^)|+m2|g°:^.(^)|)ci^<oo. 



G3 

re(0,l] • 

Note that GO means a symmetry relation Euo{x)vQ{y) = Euo{y)vo{x) that holds for an 
isotropic measure where the CFs depend only on \x — y\. Examples of such family will be 
provided later. 

Properties G0-G3 imply conditions S0-S3 for the initial measures /xo,r- Therefore, Theorem 

B implies the convergence Ht^r — ^ fico,r, t — > oo, of type (2.4). The following proposition 
means that the limiting measure /Xoo,r is close to CGD gx on the Sobolev space of distributions 
"H^'" with s,a< -n/2. 
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Proposition 4.2 Let Conditions G0-G3 hold. Then corresponding limiting measures jjioo,r 
are concentrated on any space Ti.^'°' with s,a < —n/2 and weakly converge to CCD qt on the 
space W''^: 

A*oo,r- — ^ 9t, r ^ 0. (4.6) 

Proof The convergence follows by the same arguments as in [37] from two facts (cf. Proposi- 
tions 3.2, 3.3): for any s,a with s <s < —n/2 and a <a < —n/2, 

(I) sup / \\\Y\\\l^^^^r{dY) < oo. 

re (0,1] J 

(II) For ^eV, QooA"^, ^) ^ Gri"^, ^), r ^ 0, 

where Qoo,r is the quadratic form with the integral kernel (^q^^{x — y)^, and St corresponds to 
{prix-y)^. It is important that the embedding 7i*'" C 7i*'" is compact. Property (I) can be 
checked with the help of the formula (13.3) and by using the Parseval identity: 

where fi{z) = J e-'^'{kf'dk and f2{z) = j e-'^' dk. More precisely, 

Property (I) follows from G3 as both functions fj (x) are bounded and continuous for s < —n/2. 
Furthermore, GO and (2.11) imply that = = 0, hence 

Qoo,.(*,*) = / (g^°,(x-y)*°(x),*°(t/))cixciy + I {q'J,Ax-y)<lf'{x),<lf\y))dxdy. (4.7) 
Gl and G2 together imply that 

q'^Ax-y)^T5{x-y), r ^ 0. 

Then (2.11) implies 



q 



,00 
oo.r 



Therefore, Property (II) follows from (4.7): the justification follows easily in the FT space. 
The convergence of the covariance (II) provides the convergence of the measures (4.6) as all 
the measures are Gaussian. □ 

Example Consider an initial mesure /io constructed in Example in Section 2.6.1. It satisfies 
Assumptions S0-S3 and GO. Furthermore, qQ^{z) = qQ^{z) =0, |2;| > 1 if we choose ro = 1. 
Denote by Yq{x) = {Yq{x),Yq{x)) a random function with distribution fi^. Denote by /io,r; 
r > 0, the distribution of the random function YQ^r{x) = {r^^''YQ{r~^x),r~''YQ{r^^x)) where 
ly = n/2. The corresponding CPs are gg,.(^) = r'^-^-'-^qi^^r'^z). Then all Conditions G0-G3 

hold with T := - J tTqQ^{z)dz. 
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5 Compactness of the family of measures /it 

This section gives the proof of bound (3.5). Proposition 3.2 will follow then with the help of 
the Prokhorov Theorem [37, Lemma II. 3.1] as in the proof of [37, Thm XII. 5. 2]. It is important 
that the embedding 7i C 7i~^ is compact, by virtue of Sobolev's Theorem, if £ > 0. Set: 

et = E{\u{x,t)\'^ + \Vu{x,t)\'^ + m'^\u{x,t)\'^'), xeWT. (5.1) 

The CPs of measure fit are translation invariant due to condition SI. Hence, taking expectation 
in (3.4), we get by S2, 

etl-Bfll < eo\BR+t\ < oo. (5.2) 

Here Br is the ball \x\ < R in IR"^, and \Bji\ is its volume. Taking — > oo we derive from 
(5.2) that Ct < cq: in fact, the reversibility implies then ej = cq (the mean energy density 
conservation). Hence, taking expectation in (1.5), we get (3.5): 

E\\Uo{t)Yo\\l = eo\Bn\<oo. □ 

Corollciry 5.1 Bound (3.5) implies the convergence of the integrals in (2.6). 



6 Convergence of the covariance functions 

In this section we check the convergence of the CPs of measures with the help of the PT. 
This convergence is used in Section 8. 



6.1 Mixing in terms of the spectral density 

The next proposition gives the mixing property in terms of the PT of the initial CPs . 
Assumption S2 imphes that q^ [z] is a measurable bounded function. Therefore, it belongs to 
the Schwartz space of tempered distributions as well as its PT. 

Proposition 6.1 Let the assumptions of Theorem B hold. Then q^ G L^(1R") ® M"^ , Vi, j. 
Proof Step 1 Pirst, let us prove that 

d^i^{z) e L^'(IR") ® M^, p > 1, < 2 - i - j. (6.1) 
Conditions SO, S2 and S3 imply, by [17, Lemma 17.2.3] (see Lemma 8.2 i) below), that 

\d'^qi^{z)\<Ceoip'/'{\z\), zeJR^. (6.2) 
Mixing coefficient ip is bounded, hence (6.2) and (2.10) imply (6.1): 

J \d''qo\z)\Pdz < Cel j ip'''^{\z\)dz < C^el r''-^cp^/^{r)dr < oo. (6.3) 

Step 2 By Bohner's theorem, qo = {qQ{k))dk is a complex positive-definite matrix- valued 
measure on IR"^, and S2 imphes that the total measure go(IR") is finite. On the other hand, 
(6.1) with p = 2 implies that q^^ e L^{W) (g) M^. □ 
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6.2 Proof of convergence of covariance functions 

Formulas (12.3), (12.2) and Proposition 6.1 imply for example, 

qfix -y):= e(u(x, t)»u(y,t)) (6.4) 

where the integral converges and defines a continuous function determined for all x,y E H". 
Similar integrals give a convenient modification for all functions qt'{x — y), which we will work 
with. 

Proposition 6.2 Covariance functions qt{z), i,j = 0, 1, converge for all z e IR"; 

qn^)^Qli^), t^oo, (6.5) 
where functions q^J^{z) are defined in (2.11). 
Proof. (6.4) and Proposition 6.1 imply, 

r (^) - ^(C(^) + V * q^{z)) = q'^iz), t ^ oo, (6.6) 

as the oscillatory integrals tend to zero by the Lebesgue-Riemann Lemma. For other i,j the 
proof is similar. □ 

Note that V{z) G L^(IR'^). Therefore, (6.1) with p — 1 and explicit formulas (2.11) imply 
the following 

Corollary 6.3 Functions q^ belong to L^(1R") (g) ij = 0, 1. 

Remark 6.4 A similar argument in the FT representation implies compactness in Proposition 
3.2. We provided an independent proof of the compactness in Section 5 to show the relation 
with energy conservation. 



7 Bernstein's argument for the Klein— Gordon equation 

In this and the subsequent section we develop a version of Bernstein's 'room-corridor' method. 
We use the standard integral representation for solutions, divide the domain of integration 
into 'rooms' and 'corridors' and evaluate their contribution. As a result, the value {Uo{t)YQ, \E') 
for ^ G I? is represented as the sum of weakly dependent random variables. We evaluate the 
variances of these random variables which will be important in next section. 

First, we evaluate {Y{t),^) in (3.6) by using duality arguments. For t e IR, introduce a 
'formal adjoint' operators C/o(i), U'{t) from space I> to a suitable space of distributions. For 
example, 

{Y,U[,m) = {Uoit)Y,^), ^ev, Yen. (7.1) 
Denote $(-,t) = t/o(^)^- Then (7.1) can be rewriten as 

{Y{t),^)^{Yo,^;t)), tern. (7.2) 

The adjoint groups admit a convenient description. Lemma 7.1 below displays that the action 

of groups (/^{t), U'{t) coincides, respectively, with the action of Uo{t), U{t), up to the order of 
the components. In particular, Ugit), U{t) are continuous groups of operators V V. 
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Lemma 7.1 For ^ = G V 

C/^(i)* = C/'(i)* = (7.3) 

where (f){x,t) is the solution of Eq(3.1) with the initial date {uo,vo) = (^^,^°) and ipixjt) is 
the solution of Eqn (1.3) with the initial date {uqjVq) — (^^,^°). 

Proof Differentiating (7.1) in t with Y, ^ e P, we obtain 

{Y,U^m) = {Uo{t)Y,^). (7.4) 
Group Uo{t) has the generator (3.3). The generator of U^it) is the conjugate operator 

Hence, Eqn (7.3) holds with ip = Agip. For the group U'{t) the proof is similar. □ 

Next we introduce a 'room-corridor' partition of IR". Given t > 0, choose d = dt > 1 and 
p = Pt > 0. Asymptotical relations between t, dt and pt are specified below. Set h — d + p and 

= jh, V = + d, j e Z. (7.6) 

We call the slabs Ri = {x G IR" : < < b^} 'rooms' and C{ = {x e IR" : < x"" < a^+i} 
'corridors'. Here ) , c? is the width of a room, and p of a corridor. 

Denote by Xr the indicator of the interval [0, d] and Xc that of [rf, }i\ so that I]jez(Xr(-5 ~ 
jh) + Xc{s — jh)) = 1 for (almost all) s G IR. The following decomposition holds: 

{Yo, ^■,t)) = ^((Fo, X^.*(-, ^)) + {Yo, xm-, t))), (7.7) 

where xl '■— Xr{x"' — jh) and xi '■— Xci^"' ~ jh)- Consider random variables rl, cl, where 

ri^{Yo,xm-:t)), 4-{Yo:XiH-:t)), j G Z. (7.8) 
Then (7.7) and (7.2) imply 

{Uo{t)Yo,<b) = J2{ri + cl). (7.9) 

The series in (7.9) is indeed a finite sum. In fact, (7.5) and (12.1) imply that in the FT 
representation, ^{k,t) = A'o{k)^{k,t) and ^{k,t) = g[{k)'^{k). Therefore, 

^^"^'^^ ^ (2^ / (7-10) 
This can be rewritten as a convolution 

^■,t)^nt*^, (7.11) 

where TZt — F~^Q[. The support supp^ C -B^ with an f > 0. Then the convolution repre- 
sentation (7.11) imphes that the support of the function $ at i > is a subset of an 'infiated 
future cone' 

supp$ C {(a;,t) G IR" X 1R+ : \x\<t + r}. (7.12) 
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as TZt{x) is supported by the 'future cone' \x\ < t. The last fact follows from general formulas 
(see [13, (II.4.5.12)]), or from the Paley-Wiener Theorem (see, e.g. [13, Thm II.2.5.1]), as 
^t{k) is an entire function of A; e C" satisfying suitable bounds. Finally, (7.8) implies that 

= Ct = for jh + t < —r or jh — t > r. (7-13) 
Therefore, series (7.9) becomes a sum 

{Uo{t)Yo,^) = ^{ri + cl), N,r.L (7.14) 

-Nt 

as h >1. 

Lemma 7.2 Let n > 1, m > 0, and SO— S3 hold. The following bounds hold for t > 1: 

E\ri\'<Ci^)dt/t, E\ci\'<Ci^) pt/t, jeZ. (7.15) 

Proof We discuss the first bound in (7.15) only, the second is done in a similar way. 

Step 1 Rewrite the left hand side as the integral of CPs. Definition (7.8) and Corollary 5.1 
imply by Pubini's Theorem that 

i?IWr = (x'.(a:")x'.(l/")9o(a;-|/),$(a;,t)®$(2/,t) ). (7.16) 

The following bound holds true (cf. [29, Thm X1.17 (b)]): 

sup |$(x,t)| = C(r"/2), t oo. (7.17) 

In fact, (7.10) and (12.2) imply that $ can be written as the sum 

" (2^^ / (7.18) 

where a^{uj) is a matrix whose entries are linear functions of a; or 1/uj. Let us prove the 
asymptotics (7.17) along each ray x = vt + Xq with \v\ < 1, then it holds uniformly in x e IR" 
owing to (7.12). We have by (7.18), 

<^?{vt + xo, t) = — ^ J2 I e-^('^''^'^^*-^*^^°a±(cu)^(A;) dk. (7.19) 

This is a sum of oscillatory integrals with the phase functions 4'±{k) = kv ± uj{k). Each 
function has two stationary points, solutions to the equation v — ±'Vu{k) if \v\ < 1, and has 
no if 1^1 > 1. The phase functions are nondegenerate, i.e. 

At last, ^(/c) is smooth and decay rapidly at infinity. Therefore, ^{vt + Xo,t) = C(t~"/^) 
according to the standard method of stationary phase, [14]. 

Step 2 According to (7.12) and (7.17), Eqn (7.16) implies that 

E\rl\'<Ct-'' I xi{xnho{x-y)\\dxdy = Ct-- J xiix-)dx l\\qo{z)\\dz, (7.21) 

\x\<t+r \x\<t+r R" 

where || go (-2)11 stands for the norm of a matrix (^(^{z)^. Therefore, (7.15) follows as ||go(-)ll ^ 
L^IR") by (6.1). □ 
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8 Convergence of characteristic functionals 

In this section we complete the proof of Proposition 3.3. We use a version of the CLT developed 
by Ibragimov and Linnik. If Qoo(^, ^) = 0, Proposition 3.3 is obvious. Thus, we may assume 
that for a given ^' e 

Qoo(^,*)7^0. (8.1) 



Choose < 6 < 1 and 



Lemma 8.1 The following limit holds true: 

Nt{M + i^f") + N^{<p'^\pt) + j) - 0, t^oc. (8.3) 
Proof. Function (p{r) is nonincreasing, hence by (2.10), 

r r 

r^'ip^/^r) = nj s"-V'^'(r) ds < n J s^'-^ip^/^s) ds < Op < oo. (8.4) 



Then Eqn (8.3) follows as (8.2) and (7.14) imply that Nt ^ Int. □ 
By the triangle inequality, 

\fitm-fi^m < \Eexp{z{Uo{t)Yo,^)}-Eexp{q2/t}\ 

+1 exp{-^E,^Kr} - exp{-iQoo(*, *)}| 

1 

+ \E exp{iJ2/t} - (iW{-^J2t^yt\'^}\ 
= I. + I^ + I^, (8.5) 

Nt 

where the sum J2t stands for J2 ■ We are going to show that all summands /i, I2, h tend 

j=-Nt 

to zero as t 00. 

Step (i) Eqn (7.14) implies 

h = \Eexp{i^rl}{exp{q2^ci} - 1)1 < E,^|c^l < E^l^NH'^'- (8-6) 
^From (8.6), (7.15) and (8.3) we obtain that 

Ii<CNt{pt/ty/^ ^0, t^oo. (8.7) 
Step (ii) By the triangle inequality, 

h < ^|E,^|r||^-Qoo(*,*)|<^|Qt(*,*)-Qoo(*,*)| 

= /21 + /22 + /23, (8.8) 
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where Qt is a quadratic form with the integral kernel [Qt {^iV))- Eqn (6.5) implies that 
/21 — ^ 0. As to I22, we first have that 

/22<E^I^Ml- (8-9) 

3<l 

The next lemma is a corollary of [17, Lemma 17.2.3]. 

Lemma 8.2 Let ^ be a complex random variable measurable with respect to a-algebra clA), 
f] with respect to a-algebra cr^B), and the distance dist(,4, > r > 0. 
i) Let {E\^\^y/^ < a, {E\r]\Y/^ < b. Then 

\E^7]-E^E7]\ < Cab ip^^^r). 

a) Let 1^1 < a, \r)\ < b a.s. Then 

\E^r]-E^Er]\ < Cab ip{r). 



We apply Lemma 8.2 to deduce that I22 ^ as t ^ 00. Note that rl = {Yo{x), x^(a;")(7^t * 
^)) is measurable with respect to the cr-algebra (t{RI). The distance between the different 
rooms Rl is greater or equal to pt according to (7.6). Then (8.9) and SI, S3 imply, together 
with Lemma 8.2 i), that 

I22 < CNfcp'/\pt), (8.10) 

which goes to as t ^ 00 because of (7.15) and (8.3). Finally, it remains to check that I23 0, 
t — > 00. By the Cauchy - Schwartz inequality, 

< civ,E,i^ic^f +c(E(E,w)f/'(iv,x:,£;|c^lf/^ (8.11) 

Then (7.15), (8.9) and (8.10) imply 

E{Y./if < J2E\ri\' + 2^^^^E\rlri\ < CN^dt/t + C,N,^'/\pt) < C2 < 00. 
Now (7.15), (8.11) and (8.3) yields 

l23<CiN^Pt/t + C2Nt{pt/ty/^ ^0, t^oo. (8.12) 
So, all terms /21, -^22, hs in (8.8) tend to zero. Then (8.8) imphes that 

/2<^lEt^|r^r-Qoo(*,*)H0, t^oo. (8.13) 
Step (Hi) It remains to verify that 

73 = \EeM^/t} - exp{-^E{j2^rif}\ ^0, t ^ 00. (8.14) 
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Using Lemma 8.2, ii) yields: 

\E exp{iY,/i} - n Eexp{irl}\ 

-Nt 

Nt . 

< \E exp{ir^'^^} exp{i ^ rl} — E exp{ir^'^*}E cxp{i ^ rl}\ 

-Nt+l -Nt+1 

Nt Nt 

+ \E exp{irf^^}E exp{i ^ } — ]^ ii^ exp{irf }| 

-Nt+l -Nt 

Nt Nt 

< C(p{pt) + \Eexp{t ^ ri}- [] Eexp{tri}\. 

-Nt+l -Nt+l 

We then apply Lemma 8.2, ii) recursively and get, according to Lemma 8.1, 

Nt 

\Eexp{iJ2rl} - n Eexp{iri}\ < CNMpt) ^0, t ^ oo. (8.15) 

-Nt 

It remains to check that 

Nt 2 

\l[Eexp{iri}-exp{--J2 E\rl\^}\^Q, t ^ oo. (8.16) 

-AT* 2 

According to the standard statement of the CLT (see, e.g. [24, Thm 4.7]), it suffices to verify 
the Lindeberg condition: \/e > 

Here at = Z^t-E-ktPj aiid Esf = EXsf, where Xs is the indicator of the event |/| > 6'^. Note 
that (8.13) and (8.1) imply that 

(7t ^ Qoo(*, *) 7^ 0, t^oo. 

Hence it remains to verify that Vs > 

J2E,\ri\^ 0, t^oo. (8.18) 

We check Eqn (8.18) in Section 9. This will complete the proof of Proposition 3.3. □ 

9 The Lindeberg condition 

The proof of (8.18) can be reduced to the case when for some A > we have, almost sure that 

\uo{x)\ + \vo{x)\<A<oo, xelR". (9.1) 

Then the proof of (8.18) is reduced to the convergence 

j:E\ri\'^0, t^oo (9.2) 

by using Chebyshev's inequality. The general case can be covered by standard cutoff arguments 
by taking into account that bound (7.15) for E\ri\^ depends only on cq and (p. The last fact 
is obvious from (7.21) and (6.3) with p — 1 and 7 = 0. 
We deduce (9.2) from 
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Theorem 9.1 Let the conditions of Theorem B hold and assume that (9.1) is fulfilled. Then 
for any ^ e D there exists a constant C{^) such that 

E\ri\^ <C{^)A''d^Jt\ t>l. (9.3) 

Proof. Step 1 Given four points xi, X2, X3, e IR", set: 

M^^\x^, 0:4) = E {Yo{xi) ® ... ® Yo{xi)) . 

Then, similarly to (7.16), Eqns (9.1) and (7.8) imply by the Fubini Theorem that 

i?|r^r = {k^A^i) ■ ..xiK)M!,'\xi, . . .,x,), $(xi,t) ® . . . ® <l>(x4,t)). (9.4) 

Let us analyse the domain of the integration (IR")^ in the RHS of (9.4). We partition (IR**)^ 
into three parts W2, and W^: 

4 

(IR")^ = U Wi, Wi = {x= {xi,X2, X3, Xi) e {WC'Y ■.\xi-Xi\= max \xi - Xp\}. (9.5) 

i=2 P-2,3,4 

Furthermore, given x = {xi,X2,X3,X4) e Wi, divide IR" into three parts Sj, j = 1,2,3: 
IR" = 5"! U 5*2 U 5'3, by two hyperplanes orthogonal to the segment [xi, Xi] and partitioning it 
into three equal segments, where Xi G 5*1 and Xi G 5*3. Denote by x^, Xg the two remaining 
points with p,q l.i. Set: Ai = {x E Wi : Xp E Si,Xg G S3}, Bi = {x E Wi : Xp,Xg ^ Si} 
and Ci = {x E Wi : Xp,Xg ^ S3}, i = 2,3,4. Then Wi = AiU BiU Cj. Define the function 
mo'''*(a;), x G (IR")"^, in the following way: 

(4). . j M^'^\x) -qo{xi-Xp)(^qo{xi-Xg), x e Ai, . . 

This determines mQ'^'*(a;) correctly for almost all quadruples x. Note that 

ixiixi) . . . xi{x2)qo{xi - Xp) (8) qo{xi - Xg), <^{xi,t) ® . . . ® $(x4, t)) 
= {xi{xi)xi{Xp)qoixi - Xp),^xi,t) ^ ^Xp,t)) {xi{x1)xi{x'^)qoixi - Xg),<^>{xi,t) ^ <^>{xg,t)) . 

Each factor here is bounded by C(^) dt/t. Similarly to (7.15), this can be deduced from an 
expression of type (7.16) for the factors. Therefore, the proof of (9.3) reduces to the proof of 
the bound 

It := \ixiix1) . ..xiix2Wo\^i, ■ ■ .,x,),^xi,t) . . . ® $(x4,t))| < C(*)A^d,V^^ t > 1. 

(9.7) 

Step 2 Similarly to (7.21), Eqn (7.17) implies, 

/t<C(*)t-2- j ^Xxl)...^X^l)\Yii^i\xi,...,x^)\dxidx2dx3dx^, (9.8) 

(Sf)4 

where Bl is the ball {x G IR** : < i + f}. Let us estimate m^'^-* using Lemma 8.2, ii). 
Lemma 9.2 For each i = 2, 3, 4 and almost all x G Wj the following hound holds 

\^^\xi, . . . , X4)| < CAV(ki - Xi\/Z). (9.9) 
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Proof. For x E Ai wc apply Lemma 8.2, ii) to (D^ ® (D^ = R,^ ® ]R*^-valued random variables 
^ = Yo(xi) ® Yo{xp) and r] = Yo{x,i) (g) Yo{xq). Then (9.1) implies the bound for almost all 
X G 

\m.^\x)\<CAM\x^-Xi\/3). (9.10) 

For X e Bj, we apply Lemma 8.2, ii) to ^ = ^0(2^1) and 77 = Yq{xp) (8) ^^(a^^) (8) 5^o(2^i)- Then SO 
implies a similar bound for almost all x & Bi, 

\mi^\x)\ = |M(')(x)-^ro(a;i)®^(yo(a;p)®yo(a;,)®ro(a;,))| <C'AV(ki-a;,|/3), (9.11) 

and the same for almost all x e Cj. □ 
Step 3 It remains to prove the following bounds for each i = 2, 3, 4: 



Vi(t) := J xii^i) ■ ■ ■Xi{^4)^i{^)v^{\xi - Xi\/'^)dxi dx2 dxs dx4 < Cd^t 



'2+2n-2 



(9.12) 



where Xi is an indicator of the set Wj. In fact, this integral does not depend on i, hence set 
i = 2 in the integrand: 





( \ 








j X2{x) dX4 


dxs 


dxi dx2 




W ) 







Now a key observation is that the inner integral in dx^ is 0{\xi — X2\"') as X2{x) = for 
1^4 — Xi| > |xi — X2I. This implies 



/ 



Vi{t) < Cr^ J xiix'l) J - X2\/3)\x, - Xaj" dx2 dx, J xi{x^) dx^. (9.14) 



\Bt 



The inner integral in dx2 is bounded as 

f f2{t+r) 

/ ip{\xi - X2\/3)\xi - xsl" dx2 < C{n) / r^''-^ip{r/3) dr 

Bj 

<Ci(n) sup rV^'(^/3) r^*^'V"-V'^'(^/3) (ir. 



re[0,2(t+r) 



(9.15) 



where the 'sup' and the last integral are bounded by (8.4) and (2.10), respectively. Therefore, 
(9.12) follows from (9.14). This completes the proof of Theorem 9.1. □ 

Proof of convergence (9.2). As dt < h t/Nt, bound (9.3) implies. 



□ 
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10 The scattering theory for infinite energy solutions 

In this section we develop a version of the scattering theory to deduce Theorem A from Theorem 
B. The main step is to estabhsh an asymptotics of type (1-20) for adjoint groups by using results 
of Vainberg [35]. 

Consider operators U'{t), U'^it) in the complex space H = L^{m!') ® H\M") (see (2.12)). 
The energy conservation for the KGE implies the following corollary: 

Corollary 10.1 There exists a constant C > such that V\E' G H: 

\Km\\H<c\mH, \\u'{t)nH<c\mH, tewi. (10.1) 

Lemma 10.3 below develops earher results [35, Thms 3,4,5]. Consider a family of finite semi- 
norms in H 

mlR)= J {\M^)\' + \^i{x)\' + \V^^{x)\')dx, R>0. 

\x\<R 

Denote by the subspace of functions from H with a support in the ball B^. 

Definition 10.2 He denotes the space Ui?>oi?(_R) endowed with the following convergence: a 
sequence ^„ converges to ^ in He iff > such that all ^„ e -f^(R), and ^„ converge to ^ 
in the norm || • \\{r)- 

Below, we speak of continuity of maps from He in the sense of sequential continuity. Given 
t >0, denote 

(t + l)-Mn-2(t + 2), n = 2. ^^^■^> 

Lemma 10.3 Let Assumptions El - E3 hold, and n >2. Then for any R,Ro>0 there exists 
a constant C = C{R, Rq) such that for ^ G i/(_R) 

\\U'{t)niR,)<Ce{t)mkR), t>0. (10.3) 

This lemma has been proved in [21] by using Conditions E1-E3 and a method developed 
in [36]. For the proof, the contour of integration in the /c-plane from [35] had to be curved 
loragithmically at infinity as in [36], but should not be chosen parallel to the real axis. 

The main result of this section is Theorem 10.4 below. Given t > 0, set 

ln-^(i + 2), n = 2. 

Theorem 10.4 Let Assumptions E1-E3 and S0-S3 hold, and n >2. Then there exist linear 
continuous operators W, r{t) : He ^ H such that for ^ G i^c 

U'{t)-^ = UQ{t)W^ + r{t)-^, t>0, (10.5) 

and the following bounds hold \/R > and G : 

\\rm\\H <C(i?)£i(t)l|^l|(K), t>0, (10.6) 
E\{Yo,r{tm' <C{R)el{t)mU: t>0. (10.7) 
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Proof. We apply the standard Cook method: see, e.g., [29, Thm XI. 4]. Fix ^' e and 
define W^, formally, as 

W-^ = lim U'{-t)U'{t)^ = ^ + / —U'Q{-t)U'{t)^ dt. 
*->oo J dt 



We have to prove the convergence of the integral in norm in space H. First, observe that 

where A!q and A! are the generators to groups U^it), U'{t), respectively. Similarly to (7.5), we 
have 

^'=(J o)' 

n 

where A = {dj — iA^Y — m?. Therefore, 

jU',{-t)U[m = U^{-t){A' - X)U'm. (10.9) 
Now (10.8) and (7.5) imply 

n 

Furthermore, E2 implies that L = J2 (dj — iAjy — A is a first order partial differential operator 
with the coefficients vanishing for \x\ > Rq. Thus, (10.1) and (10.3) imply that 



\m~t){A'-x)u'{m\H < c\\{A'-x)u'm\\H = c\\[iA'-x)u'm) 

< C, ||(t/'(t)f)'|ki(s«„) < C{R)e{t)m(R), t > 0. (10.10) 

Hence (10.9) impfies 

I \\jfl{-t)U'{t)nHdt < C{R)ei{s)m^R), s > 0. (10.11) 

s 

Therefore, (10.5) and (10.6) follow by (10.1). It remains to prove (10.7). First, similarly to 
(7.16), 

E{Yo, r(i)*)2 = {qo{x - y), r(i)*(x) r(i)*(t/) ) (10.12) 
Therefore, the Shur Lemma implies (similarly to (7.21)) 

E{Yo,rit)^)^ < llgolUi MmWv^ Wrmh^, (10.13) 
where the norms || • \\lp have an obvious meaning. Finally, (10.6) implies for ^ e 

WrmU^ < C\\rit)^\\H < CiR)eiit)\M{R), (10-14) 

Therefore, (10.7) follows from (10.13) since ||?o||li < oo by (6.1). □ 
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11 Convergence to equilibrium for variable coefficients 

The assertion of Theorem A follows from two propositions below: 

Proposition 11.1 The family of the measures {nt, t G IR}, is weakly compact in , Vs > 0. 
Proposition 11.2 For any ^ e D 

fit{^) = J cxp{t{Y, ^)) fn{dY) ^ exp{~Q^{W^, )}, t ^ oo. (11.1) 

We deduce these propositions from Propositions 3.2 and 3.3, respectively, with the help of 
Theorem 10.4. 

Proof of Proposition 11.1 Similarly to Proposition 3.2, Proposition 11.1 follows from the 
bounds 

sup E\\U{t)Yo\\R < oo, R>0. (11.2) 
For the proof, write the solution to (1.3) in the form 

u{x,t) = v{x,t) + w{x,t). (11-3) 
Here v{x, t) is the solution to (3.1), and w{x, t) is the solution to the following Cauchy problem 
' w{x, t) = ELi(^fe - iAk{x)yw{x, t) - rr? w{x, t) 

-E'U^^M^)dkv{x,t) -Ek=ii^9kAk{x) + Al{x))v{x,t), (11.4) 

w\t=o = 0, ^14=0 = 0, xe IR". 
Then (11.3) implies 

E\\U{t)Yo\\R< E\\Uo{t)Yo\\R + E\\{w{;t),w{;tmR- (11-5) 
By Proposition 3.1 we have 

sup E\\Uo{t)Yo\\R < oo. (11.6) 

t>o 

It remains to estimate the second term in the right hand side of (11.5). The Duhamel repre- 
sentation for the solution to (11.4) gives 

t 

{w,w) = J U{t- s){0,il){-,s))ds, (11.7) 



n n 

where iplxjs) = —2i J2 Ak{x)dkv{x, s) — J2 {'idkAk{x) + Al{x))v{x, s) . Assumption E2 implies 

k=l k=l 

that supp'0(-,s) C Brq. Moreover, 

||(0,V(-,«))lko < C\\v(.,s)\\h^^s^^) < C\\Uo(s)Yo\U. (11.8) 

The decay estimates of type (10.3) hold for the group U (t), as well as for U'{t), as both groups 
correspond to the same equation by Lemma 7.1. Hence, we have from (11.8), 

\\U{t - s){Q,^{; s))\\r < C{R)e{t - s)||(0, V^(-, s))||^„ < C,{R)e{t - s)\\Uo{s)Y4r,, (11.9) 
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where e{-) is defined in (10.2). Therefore, (11.7) and (11.6) imply 

t 

E\\{w{-,t),w{-,t))\\R<C{R) J e{t - s)E\\Uo{s)Yo\\R,ds < C2{R) < oo, t>0. (11.10) 



Then (11.6) and (11.5) imply (11.2). □ 
Proof of Proposition 11.2 (10.5) and (10.7) imply by Cauchy-Schwartz, 

\Eexpi{U{t)Yo,^) - Eei^pt{Yo,U;,{t)W^)\ < E\{Yo,r{t)^)\ 

It remains to prove that 

Eexpi{Yo,U^{t)W^) ^ exp{-^Qoo(W^*,W^* )}, t ^ oo. (11.11) 

This docs not follow directly from Proposition 3.3 since generally, W"^ ^ V. We approximate 
W"^ by functions from V. W"^ G H, and V is dense in H. Hence, for any e > there exists 
$ e X> such that 

\\W^-^H<e. (11-12) 
Therefore, we can derive (11-11) by the triangle inequality 

|£;expz(ro,t/oWW^*) -exp{-^Qoo(W^*,W^* )}| 
< \Eexpi{Yo,U;,(t)W^) - Eexpi{Yo,U;,{t)^)\ 

+£;|expi(C/o(t)yo,^) -exp{-^Qoo(^,^ )}| 

+1 exp{-^Qoo(*, * )} - eM — QooiW^, )}\. (11.13) 
Applying Cauchy-Schwartz, we get, similarly to (10.12)-(10.14), that 

E\{Yo, U[,{t){W^ - $))! < (E|(Fo, t/^(t)(W^* - < C\K{t){W^ - ^)\\h. 

Hence, (10.1) and (11.12) imply 

E\{Yo,U^{t){W^-^))\<Ce, t>0. (11.14) 

Now we can estimate each term in the right hand side of (11.13). The first term is 0{e) 
uniformly in i > by (11.14). The second term converges to zero as i — > oo by Proposition 
3.3 since $ e P. Finally, the third term is 0{e) owing to (11.12) and the continuity of the 
quadratic form Qoo(^5 ^) in i^^(lR") ® (D^. The continuity follows from the Shur Lemma since 
the integral kernels q^c^{z) G L^iMT') ® by Corollary 6.3. Now the convergence in (11.11) 
follows since e > is arbitrary. □ 

12 Appendix A. Fourier transform calculations 

Consider the covariance functions of the sohitions to the system (3.2). Let F : w ^ w denote 
the FT of a tempered distribution w G S"(IR") (see, e.g. [13]). We also use this notation for 
vector- and matrix-valued functions. 
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12.1 Dynamics in the FT space 

In the FT representation, the system (3.2) becomes Y{k,t) — Ao{k)Y{k,t), hence 

Y{k, t) = gt{k)Yo{k), Gtik) = exp{Ao{k)t). 

Here we denote 



/ 

Mk) = I , Gt{k) = 




/ smcut \ 
cos ujt 

y —CO sin cut cos cut j 



12.2 Covariance matrices in the FT space 
Lemma 12.1 In the sense of matrix-valued distributions 

qt{x - y) := E{Y{x,t) ® Y{y,t) ) = F^X-ySt{k)qo{k)G't{k), t G IR. 
Proof Translation invariance (1.8) imphes 

e{Yo{x) ®c Yoiv) ) = C^{x - y), e(Yo{x) ®c ) = Co"(^ - v), 

where ®c stands for tensor product of complex vectors. Therefore, 

E[%{k) 0c %{k') ) = F^^kFy^k' C+{x -y) = {2^T5{k + k')C^{k), 
E[%{k) ®cY^)) = F,^kFy-.-k'C^{x -y)^ {2T:Y5{k - k')C^{k). 

Now (12.1) and (12.2) give in the matrix notation that 

E{Y{k,t) ®cY{k',t) ) = {2T:Y5{k + k')gt{k)C^{k)g[{k), 
E[Y{k,t) ®cY{k',t) ) = (27r)"5(A; - k')gt{k)C^ {k)g',{k) . 

Therefore, by the inverse FT formula we get 

E{Y{x,t) ®cY{y,t) ) = F^l^_^g,{k)C^{k)g',{k), 

e{y{x, t) ®c nyJ) ) = F^l,_^gt{k)c^{k)g[{k). 

Then (12.3) follows by linearity. 

13 Appendix B. Measures in Sobolev's spaces 

Here we formally verify bound (4.5) for s, a < —n/2. Definition (4.2) implies for u 



ix{k—k') I 



«lli«= /(^) I e-'-^'-'KkY{kru{k)u{k')dkdk' 



dx. 
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Let fi{du) be a translation-invariant measure in with a CF Q{x,y) = q{x — y). Similarly 
to (12.5), (12.4), we get 

j u{k)^k^n{du) = (27r)"5(A; - k') trg(A;). (13.2) 

Then, integrating (13.1) with respect to the measure ii{du), we get the formula 

/ MlMdu) = / i^r^'dx J (k)'^ trq{k) dk. (13.3) 

Applying it to q{k) — T with a,s < —n/2 and to q{k) — T{k'^ + m^)~^ with 1 + s instead of 
s, we get (4.5). 
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